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Using a new rotating SENSitivity Encoding (rotating-SENSE) algorithm, we have successfully demon- 
strated that the rotating radiofrequency coil array (RRFCA) was capable of achieving a significant reduc- 
tion in scan time and a uniform image reconstruction for a homogeneous phantom at 7 Tesla. However, at 
7 Tesla the in vivo sensitivity profiles (BĮ ) become distinct at various angular positions. Therefore, sensi- 
tivity maps at other angular positions cannot be obtained by numerically rotating the acquired ones. In 
this work, a novel sensitivity estimation method for the RRFCA was developed and validated with human 
brain imaging. This method employed a library database and registration techniques to estimate coil sen- 
sitivity at an arbitrary angular position. The estimated sensitivity maps were then compared to the 
acquired sensitivity maps. The results indicate that the proposed method is capable of accurately estimat- 
ing both magnitude and phase of sensitivity at an arbitrary angular position, which enables us to employ 
the rotating-SENSE algorithm to accelerate acquisition and reconstruct image. Compared to a stationary 
coil array with the same number of coil elements, the RRFCA was able to reconstruct images with better 
quality at a high reduction factor. It is hoped that the proposed rotation-dependent sensitivity estimation 
algorithm and the acceleration ability of the RRFCA will be particularly useful for ultra high field MRI. 


© 2014 Elsevier Inc. All rights reserved. 


1. Introduction 


Reducing the scan time of magnetic resonance imaging (MRI) is 
important for both clinical diagnosis and scientific research. Clini- 
cally, acquiring high resolution images normally requires long scan 
duration, which may discomfort patients and generate motion- 
related artifacts that affect diagnostic accuracy. In research, espe- 
cially in functional MRI (fMRI), shorter scan time facilitates 
improving quantification of the brain activity. The fast full k-space 
acquisition schemes have been developed for modern MRI systems 
to shorten scan time, such as echo planer imaging (EPI) [1], fast 
spin echo [2] and fast gradient echo imaging [3]. More recently, 
phased array coils (PACs) [4,5] have facilitated partial k-space 
acquisition to further reduce the scan time. Images can then be 
reconstructed with different algorithms, such as simultaneous 
acquisition of spatial harmonics (SMASH) [6], SENSitivity Encoding 
(SENSE) [7,8], generalized autocalibrating partially parallel acquisi- 
tions (GRAPPA) [9]. 
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As an alternative method, we have introduced the rotating 
radiofrequency coil (RRFC) [10,11] and the rotating radiofrequency 
coil array (RRFCA) [12-14] to reduce scan duration by using 
SENSE-like reconstruction methods. Both rotating RF systems are 
capable of encoding k-space data with varying sensitivity profiles 
(B; ), which benefit the signal encoding process and reconstruction 
algorithms. By adopting the time-division multiplexing sensitivity 
encoding (TDM-SENSE), the RRFC demonstrated a 2-fold scan time 
reduction for human brain imaging [11]. The RRFCA used the rotat- 
ing-SENSE algorithm that featured time-varying sensitivity profiles 
for acquisition acceleration and image reconstruction. A 4-element 
RRFCA was prototyped and tested with a homogeneous phantom 
at 7 Tesla, which exhibited improved g-maps [7] and uniform sig- 
nal-to-noise ratio (SNR) distributions compared with a traditional 
stationary 8-element coil array [12,14]. 

In this work, the feasibility of using the RRFCA for human brain 
imaging was tested. Since the rotating-SENSE algorithm needs sen- 
sitivity information for image reconstruction, in our previous stud- 
ies [11,14], a small number of sensitivity maps acquired at several 
angular positions were numerically rotated to estimate sensitivity 
at other positions. This technique was used for human brain imag- 
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ing at 2 Tesla [10,11] owing to the negligible coil-tissue interaction 
[15-19] (dielectric resonance [20-22], or constructive and destruc- 
tive interference [23,24]). This approach was also used for imaging 
a homogeneous phantom at 7 Tesla, provided the phantom and the 
RRFCA system are both symmetrical [14]. However, this method is 
not applicable for in vivo sensitivity estimation at 7 Tesla, because 
the in vivo sensitivity maps at different angular positions vary sig- 
nificantly. Since measuring sensitivity at every angular position is 
impractical and compromises the fast imaging purpose, a practical 
and robust sensitivity estimation method is needed for in vivo 
applications of the RRFCA at ultra high fields. 

Here we investigate a novel sensitivity estimation method spe- 
cially designed for in vivo applications of the RRFCA. Instead of sim- 
ply rotating the acquired sensitivity to new angular positions, the 
in vivo sensitivity maps will be deformed in a non-linear fashion. 
These deformations can be calculated with image registration tech- 
niques by registering in vivo sensitivity to sensitivity in the library 
which was acquired from scans of volunteers or numerical calcula- 
tions. This approach is based on the observation that the B; map is 
not particularly sensitive to small local changes, in terms of dielec- 
tric properties and structures of tissues. Instead, B; maps are typ- 
ically related to the global dielectric property distribution relative 
to the RF system [24-26], with lower spatial frequencies. This sen- 
sitivity deformation can be modelled numerically using image reg- 
istration techniques [27-29], as these techniques are commonly 
used in finding the spatial correspondence between two images. 
Importantly, image registration has recently found application in 
modelling the magnetic field variations due to changes in dielectric 
distributions between different subjects [30]. 

In this work, the proposed numerical method will be applied to 
estimate the rotated sensitivity maps, which will then be used to 
optimise the rotating technique and reconstruct the images. The 
optimisation of the rotating scheme is achieved by minimising 
the maximum g-factor, which suppresses the noise amplification 
in the imaging acceleration. The g-maps, reconstructed images 
and SNR maps of the RRFCA and stationary coil array will then 
be compared. 


2. Methods and materials 


As described in previous work [14], the RRFCA moves to various 
angular positions during acquisition, so that different k-space 
phase-encoding lines are associated with distinct coil sensitivity 
profiles which improve encoding capability. In order to reduce 
the overall scan time, sensitivity maps at most angular positions 
are estimated from a small number of measured in vivo sensitivity 
maps, by employing the proposed algorithm. The linear superposi- 
tion of individual sensitivity maps with global coverage of the sam- 
ple can benefit the registration algorithm [31] (details are 
discussed in later sections). However, four physical elements with 
90° intervals of the RRFCA prototype are insufficient to provide a 
complete coverage. Consequently, two sets of sensitivity maps, 
with angular separation of 45° (e.g. position 1: 0° + 45° in Fig. 1), 
are used together to simulate a rotating array with 8 elements 
equidistantly distributed in angular direction. Considering that 
all elements are identical, the range of angular displacement that 
the RRFCA needs to travel is from —22.5° to +22.5°. In order to test 
the robustness of the estimation algorithm at the maximum angu- 
lar rotation, sensitivity with 22° displacement (position 3) from 
position 1 is estimated. In addition, the intermediate angular posi- 
tion displacement (10° at position 2) is also estimated. (as shown 
in Fig. 1, estimate B} at positions 2 and 3 from acquired B} at posi- 
tion 1). To verify the robustness of the proposed algorithm, two 
volunteers of different genders and distinctly different head sizes 
were imaged. 


position] = 0°+45° 
position2 = 10°4+55° 
position3 = 22°4+67° 


Fig. 1. The sensitivity profiles of the 4-element RRFCA were acquired at angular 
positions at 0°, 10°, 22°, 45°, 55° and 67°, denoted with red, yellow, green, dashed 
red, dashed yellow and dashed green. Sensitivity maps acquired at these positions 
were recombined as three positions for better registration performance: position 1 
(0° + 45°), position 2 (10° + 55°), position 3 (22° + 67°). (For interpretation of the 
references to color in this figure legend, the reader is referred to the web version of 
this article.) 


2.1, Registration based in vivo sensitivity estimation 


In this work, the new sensitivity estimation method uses the 
acquired sensitivity at the initial position to estimate the sensitiv- 
ity at other positions by employing a sensitivity library and regis- 
tration techniques. The library provides source images, which are 
made up of sensitivity maps at all angular positions acquired dur- 
ing scans of volunteers. The registration techniques are employed 
to find the spatial transformation that aligns the source image 
(library sensitivity) to the target image (actual acquired sensitivity) 
from its initial position, and this transformation is applied to other 
angular positions for estimating corresponding sensitivity maps. 
The rotational sensitivity at arbitrary angular positions is acquired 
by following four steps: 


(1) Create a library by acquiring sensitivity maps at all angular 
positions from the scans of the volunteers. During the actual 
patient scanning, sensitivity maps at an initial position (i.e., 
position 1 in Fig. 1) are acquired. Image processing is applied 
to both groups of sensitivity maps to smooth profiles and 
correct singular values. 

(2) Instead of registering individual coil sensitivity maps, they 
are linearly combined before being registered, to improve 
accuracy and efficiency. The combination coefficients are 
determined with a condition number of 1 to achieve the 
optimal sensitivity estimation when individual sensitivity 
profiles are later extracted (see Step (4)). a 

(3) As illustrated using a flow chart in Fig. 2, the source profile S$ 
(combined 8 library B} maps at position 1) is registered to 


the target profile TÌ (combined 8 in vivo B] maps at position 


1). The corresponding transformation ¥% is extracted and 
applied to the combined library sensitivity map at arbitrary 
angular position («°) to obtain the estimated combined sen- 


bats: = 
sitivity map E,. 
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Repeat j times, j = 8 


Fig. 2. Flow chart of the registration based rotation-dependent sensitivity estimation. 


(4) Repeat step (3) 8 times for all the linear combinations to 
estimate all the combined sensitivity maps. The individual 


sensitivity map E} at «° angular position (i.e., positions 2, 3 
in Fig. 1) is calculated by multiplying the inverse of the coef- 
ficient matrix. N? in Fig. 2 denotes the number of pixels. 


2.1.1. In vivo sensitivity mapping and singular value correction 

The most common sensitivity mapping method for a stationary 
coil array is to derive relative coil sensitivity by dividing the indi- 
vidual coil image using a predominantly uniform reference image. 
Both the coil image and the reference image can be obtained from 
either full k-space sampling [7,32] or fully sampled central k-space 
[33,34]. Without a uniform volume transmit coil, the reference 
image is typically approximated as the root sum of square (RSS) 
image: 


J +2 
mss =| ya (1a) 
ae 
Si = Rss es 


where J is the total number of coil elements in an array; Ij is the full- 
FOV image obtained with the j-th coil; S; is the sensitivity map of 
the individual coil. 

The direct division in Eq. (1b) can easily generate singular val- 
ues. The different singular value distributions of each coil can cause 
inaccurate reconstructions which, particularly in this work, can 
misinform the registration process. 

Singular values are commonly seen at the interfaces (i.e. skin 
and skull, skin and air) for both magnitude and phase of the sensi- 
tivity maps. In addition, the phase of sensitivity is often unpredict- 
ably wrapped for the ultra high field MRI. To prepare sensitivity 
profiles for registration, a multi-level singular value removal algo- 
rithm was developed. Similar to methods employed in [35], the 
map was first divided into reliable and unreliable regions. The 
unreliable regions consist of singular values of the magnitude 
map, sensitivity voids and phase discontinuities [36]. An interpola- 
tion/extrapolation procedure [35], based on polynomial fitting, was 
then performed to correct the sensitivity profiles in the unreliable 
regions. 


2.1.2. Optimal sensitivity combination 

Due to the complex coil-tissue interactions at ultra high fields, 
after loading a heterogeneous subject like the human brain, the 
sensitivity of each coil becomes distinct [25]. For this reason, single 


Estimated 
Sensitivity 
EL 

(8xN?) 


coil sensitivity based estimation is compromised by insufficient 
coil-tissue information. In addition, the registration algorithm 
works less efficiently with signal voids that are commonly seen 
in single-channel sensitivity maps. For better estimation, sensitiv- 
ity profiles are linearly combined before applying registration. A 
similar method has been reported for optimal B,; mapping [31]. 
This process is expressed as: 


Ccoef Ssingle = Scombined (2) 


where Ccoer denotes an 8 x 8 coefficient matrix; Ssingie is the sensi- 
tivity matrix resized to 8 x N? (resolution is N x N) and each row 
is a single coil sensitivity. Scombinea Contains the linearly combined 
sensitivity maps. 

During rotating, the coil-tissue interaction is fully transformed 
from one angular position to another and Eq. (2) is applied accord- 
ingly. As shown in Fig. 2, to decompose the estimated individual 


coil sensitivity E! from estimated combined sensitivity maps E$, 
the inversion of Ceoef is multiplied to both sides of equation. With 
the method described in [31], the matrix Ceoef is determined with 
condition number 1 to minimise error. 


2.1.3. The registration technique for sensitivity estimation 

The goal of registration is to find the spatial correspondence 
between the source and target images. Registration techniques 
are widely used in medical imaging [27,28,37,38], such as when 
information from multiple modalities are combined (computed 
tomography (CT), MRI, positron emission tomography (PET)), 
intra-subject motion correction, distortion correction, dynamic 
imaging reconstruction and high-field MRI safety assessment 
[30]. A typical registration algorithm includes two components: a 
similarity metric and geometric deformation. The registration algo- 
rithm employed here followed, in part, previously published works 
[29,39], while various components were adapted to suit the cur- 
rent application. For example, a similarity metric was modified to 
cope with complex-valued sensitivity maps. 


2.1.3.1. Similarity metric. The similarity metric measures how well 
the source image aligns to the target image. Various metrics have 
been developed for image registration, such as squared difference 
(SD), mutual information (MI) [40,41] and pattern intensity (PI) 
[42]. 

Sensitivity maps are complex images. Both magnitude and 
phase of By and By exhibit a rotational property [25] which offers 
opportunities to apply image registration techniques for estimat- 
ing sensitivity maps. Magnitude images are relatively easy to be 
registered because the features are clear; however, registering 
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phase images is problematic due to phase-wrapping. Even when 
the phase images are unwrapped, the unwrapping quality and 
phase ranges are difficult to control which can undermine the reg- 
istration performance. To avoid this problem, the registration pro- 
cess in this work was not solely performed on phase images, but on 
the complex sensitivity maps. The PI metric was modified to cope 
with complex numbers. The SD metric was chosen for registering 
magnitude images due to its fast processing speed. Two registra- 
tion processes were set to different convergence parameters for 
optimal performance. 

For a source image I,(x,y) and a target image I[;(x,y), a two- 
dimensional SD and PI are defined as: 


C Ea Dya y) = he(x,y))? 


SD XY 


(3a) 


o2 
a? + (laip (XY)-laip (u,v)? 
—v)?<r2 


P-D 
XY (x—u) +(y: 


laip = 1 — Is 


(3b) 


where x and y are the voxel coordinates in the image, X and Y are 
the numbers of the voxels along each dimension. In Eq. (3b), [air 
is the difference image. The parameter r defines the size of the 
neighbourhood, in which the variations are taken into account. 
The parameter ø is a sensitivity controller to decide whether the 
variation is a structure or not. 


2.1.3.2. Geometric deformation. Two types of geometric deforma- 
tion, affine transformation and non-rigid transformation, were 
employed in this work. The affine transformation included transla- 
tions and rotations, making it best suited for global transforma- 
tions. The non-rigid transformation was applied to compensate 
the inefficiency of the affine transformation in local areas. In the 
non-rigid transformation, a deformation field recorded all dis- 
placement vectors at each voxel from an aligned source image to 
the target image. In this work, a free-from deformation (FFD) based 
on B-splines [43] was adopted for non-rigid transformations. The 
FFD worked by deforming a source image by manipulating a mesh 
of control grids [44]. To improve its efficiency, a coarse-to-fine 
multilevel B-spline approximation was adopted to generate a ser- 
ies of B-spline functions incorporating bicubic interpolation func- 
tions for calculation of the deformed pixel values. 


transmit array @ 400mm 


RRFCA Ø 280mm 


(a) 


(b) 
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2.2. Image reconstruction and encoding optimisation 


The rotating coil array emulates a large amount of coils with 
only four physical RF elements, thereby improving the condition 
of the encoding matrix to facilitate a higher reduction factor. How- 
ever, the encoding ability of RRFCA can be further improved by 
strategically choosing the rotating degree for each stepping, which 
is determined by the g-map based optimisation algorithm [14]. 


2.2.1. Image reconstruction with rotating-SENSE 

The number of sensitivity profiles limits the sensitivity encod- 
ing ability of a stationary coil array; therefore, the sensitivity pro- 
files in the encoding matrix remain unchanged from row to row for 
each coil. However, the encoding matrix (Ag) of the RRFCA has 
more variations in its rows by taking rotational sensitivity into 
account: 


ES, 


ES) 
Ar = 


ES, 
where E = e"? and S; = S;(rp, t). 


K € [1, M x N/R], t € [1, MIR], N is the length of a k-space line, M is 
the number of k-space lines and R is the reduction factor. E and S; 
denote the Fourier encoding and sensitivity encoding matrices, 
respectively. K, denotes the x-th k-space sampling position; rp 
and t in S; denote the sensitivity at the position of p-th voxel of 
the step t. 

The Fourier encoding kernel is consistent between the station- 
ary and rotating array; however, sensitivity maps for individual 
coils j at each step t are different as the coil rotates. This variation 
improves the condition of encoding matrix and can be exploited for 
further scan time reduction [14]. 

The noise behaviour analysis of the rotating-SENSE is similar to 
that of the traditional stationary array. However, to employ g-maps 
[7] for noise analysis, the sensitivity matrix S needs to be delin- 


eated from encoding matrix Ag: 
Ar = FS (5) 


where F and S represent the Fourier encoding matrix and sensitivity 
matrix, respectively. The sensitivity matrix S consists of the 


Fig. 3. (a) Top view of the RRFCA setup. The outside layer is the stationary 8-element transmit coil array with 400 mm diameter, inner layer is the 4-element RRFCA with 
280 mm diameter. (b) RRFCA system loaded with a volunteer. [A] RRFCA system. [B] Volunteer. [C] Extended shaft pointing to rear of MRI. (c) Close-up of the RRFCA system. 
[D] Transmit coils made of copper patches attached to outside stationary former. [E] Receive coil attached to the internal rotatable former. [F] Detuning circuits on both 


transmit and receive coils. 
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acquired sensitivity and estimated sensitivity profiles. See [14] for 
the process of separating F and S. 


2.2.2. Sensitivity encoding optimisation 

The proposed method of estimating in vivo sensitivity maps 
enables us to investigate the optimal rotating scheme using numer- 
ical simulations. By optimally choosing the sampling position of 
each step, the sensitivity encoding capability can be maximised. 
At each angular position, the coil sampled one phase-encoding line. 
The angular displacement 0 between adjacent acquisitions was 
determined by achieving the best imaging acceleration perfor- 
mance as follows: 


arg nin Er Dy- Bry(6) - 2 


MxN 
where 0 denotes the angular displacement between two k-space 
lines, Zx y(0) is the g-factor calculation [7] at voxel (x, y) in an image 
with M x N pixels. 


(6) 


2.3. Experimental validation 


The 4-channel RRFCA was used to scan two subjects, A and B, at 
6 angular positions (0°, 10°, 22°, 45°, 55° and 67° in Fig. 1). These 
data were firstly combined to simulate an 8-channel RRFCA (for 
reasons described previously). To test the proposed algorithm in 
an intra-subject case, images of subject A were acquired at two 
slices with 12 mm spacing. Sensitivity maps of one slice were used 
as the library. In the inter-subject case, sensitivity maps of subject 
A were used as the library to estimate sensitivity maps of subject B. 
The image registration techniques were then applied to both data- 
sets for the estimation of the rotation-dependent sensitivity pro- 
files at desired positions. 

To study the rotating coil array in the reception mode indepen- 
dently without the complication from changing transmission pro- 
files, a 4-element RRFCA [14] was used as a receive-only coil array. 
An additional transmit coil array was built to provide an 
unchanged and relatively uniform transmission. The experimental 
setup is shown in Fig. 3(b) with one healthy volunteer on a whole- 


body 7T MRI scanner (Magnetom 7T MRI, Siemens Medical Solu- 
tions, Erlangen, Germany). All imaging protocols were approved 
by institutional review board of Institute of Biophysics of Chinese 
Academy of Sciences (Beijing), and signed consent forms were col- 
lected. Small-flip-angle (FA = 30°) GRE images of two slices were 
acquired with TE/TR=4/1000ms (in plane voxel size= 
2mm x 2 mm). 

The RRFCA system consisted of three layers. The stationary 
inner layer was isolated from the rotating coils and provided sup- 
port to a patient’s head. The receive coils were attached to a rotat- 
able middle layer and the transmit coils were attached to a 
stationary outer layer. The radius of the inner layer was set to 
125 mm to comfortably accommodate the human head. As shown 
in Fig. 3(a) and (c), the transmit coil array consisted of eight loop 
coils, each of which was 160 mm in length and 130 mm in width. 
They were attached equidistantly to a coil former. A capacitance 
decoupling method was used to keep the coupling at a reasonably 
low level of —12 dB ~ —14 dB when loaded. To reduce the interac- 
tion and parasitic capacitance between the transmit and receive 
coils, besides employing active detuning circuits, a transmit coil 
former with a large radius should be used. However, in order to 
minimise the signal drop caused by the increased distance and to 
guarantee consistent rotation, the radius of the middle layer was 
set to 140 mm to provide a 15 mm separation from the inner layer. 
(Fig. 3(a), Ø transmit? 400 mm, Ø receive! 280 mm). 

An extended shaft [14] was used to adjust the rotation angle 
outside of the MRI tunnel without repositioning the coil-subject 
set for acquisition at each angular (stepping) position. However, 
unlike the experiment setup using phantoms [14], the rotation 
indicator and the extended shaft were installed at the rear of the 
magnet bore (Fig. 3(b)) to allow enough space for the patient. 


3. Results 


Fig. 4(a) and (b) shows the gradient recalled echo (GRE) images 
of two slices from subject A. The GRE image of subject B (same slice 
location as Fig. 4(b)) is shown in Fig. 4(c). Their corresponding sen- 
sitivity maps are shown in (d)-(f). Since there is a large difference 
in head size and thus in global dielectric properties between 


Fig. 4. (a) and (b) are GRE images of two slices (12 mm separation) from subject A at position 1 (0°). (c) Image from subject B of the same slice position in (b) at position 1 (0°). 


Their corresponding sensitivity maps are shown in (d)-(f) respectively. 
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Fig. 5. (a) and (c) are raw magnitude and phase plots of sensitivity. The areas in the red boxes are enlarged to show singular values. In (b) and (d), the red dots denote the 
singular values. Compared to adjacent areas, these spikes have higher gradients and are mostly seen at the interfaces. Besides singular values, undulating errors across the 
raw sensitivity profiles also need to be corrected for better registration efficiency. (For interpretation of the references to color in this figure legend, the reader is referred to 


the web version of this article.) 


subjects, large sensitivity variations are found even at the same 
slice location. As seen from Fig. 4(a) and (b), the anatomical struc- 
ture and dielectric property of the two slices (subject A) from one 
subject are very different, but their sensitivity maps have minor 
changes shown in Fig. 4(d) and (e). This correlates with the obser- 
vation that the sensitivity is more related to global dielectric 
changes rather than local changes. 

The raw magnitude and phase maps of sensitivity derived from 
Eq. (1a,b) are shown in Fig. 5(a) and (c), respectively. As predicted, 
coil sensitivity maps should be smooth with a small local gradient. 
However, in Fig. 5(b) and (d), we can see that the coil sensitivity 
maps had abrupt changes which were associated with high gradi- 
ents (marked as red dots). Magnified regions in the red boxes show 
the singular values. In addition, the phase map in Fig. 5(c) also 
exhibits the phase wrapping at the skin/air and the skull/tissue 
interface. 

With the developed multi-level fitting algorithm, both the mag- 
nitude and phase images were smooth with the singular values 
corrected. In Fig. 6(a), the signal voids were also extrapolated for 
image registration purposes. Compared to Fig. 5(b) and (d), the 
high gradient and undulating errors were corrected. Both the mag- 
nitude and phase plots of sensitivity were smooth and natural. 


The experimentally acquired and numerically estimated sensi- 
tivity maps for the intra-subject case are compared in Fig. 7. With 
a 10° angular displacement (p2 in Fig. 7(a)), the estimated magni- 
tude maps in the second row of Fig. 7(a) were very similar to the 
measured sensitivity maps in the first row. Both the global features 
and the local details were estimated accurately. Estimated phase 
maps in the fourth row only showed small local variations com- 
pared to the acquired maps. Combining the magnitude and phase 
plots into a complex-numbered sensitivity, the root-mean-square 
error (RMSE) of the estimated sensitivity is 0.048. In Fig. 7(b), the 
angular displacement increased to 22° (position 3), which was also 
the largest displacement for the rotational sensitivity estimation 
under the current configuration. From the comparisons between 
the estimated and measured sensitivity maps in Fig. 7(b), we note 
that the minor local discrepancies started to present, although the 
global features of sensitivity maps were well captured by the reg- 
istration based algorithm. The RMSE increased to 0.068 at p3, 
which suggested a slightly decreased accuracy of estimating sensi- 
tivity with a larger angular displacement. 

For the inter-subject case, the experimentally acquired and 
numerically estimated sensitivity maps at position 2 (p2 in 
Fig. 8(a)) and position 3 (p3 in Fig. 8(b)) are compared in Fig. 8. 
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Fig. 6. (a) magnitude (top row) and phase (bottom row) plots of refined sensitivity map at 0°. (b) 3D magnitude (top row) and phase (bottom row) plots of refined sensitivity 


map for single coil. 


The female subject had a noticeably smaller head than that of the 
male subject, therefore the sensitivity maps of two subjects have 
more differences compared to the intra-subject case. However, as 
shown in Fig. 8(a), the registration based algorithm was able to esti- 
mate the magnitude and phase maps of sensitivity accurately, both 
globally and locally. The RMSE was 0.054, which is slightly higher 
than that of the same position in the intra-subject case 
(RMSE = 0.048). Different from the intra-subject case, the RMSE 
reduced to 0.045 when the angle was then increased to 22° for posi- 
tion 3. 

To illustrate the imaging acceleration and reconstruction ability 
of the RRFCA, the g-maps, reconstructed images and the SNR maps 
are shown in Fig. 9 when the reduction factor was four (R = 4). The 
measurements with 4-element RRFCA sampled at two angular 
positions (position 1: 0° + 45°) were used together to emulate an 
8-element stationary coil array. In the first and second rows, 
results reconstructed from the experimental data for 4- and 8-ele- 
ment stationary arrays are shown, respectively. However, the coil 
elements of the RRFCA were naturally decoupled (—18 dB S,,), 
which would be more difficult to realise with an 8-element coil 
array in the same size. Therefore the max g-factor in the second 
row may be higher than 1.7 and the averaged relative-SNR may 
be lower than 2.35 in practice. 


Since the SNR is hard to calculate accurately for ultra high field 
MRI, the SNR maps in this work were calculated relative to that of a 
single rotating coil without undersampling (R = 1). As we can see in 
Fig. 9, with a high reduction factor (R = 4), the 4-element stationary 
coil array had a very high g-factor (max-g = 3.7) and a 0.95 aver- 
aged relative-SNR, which corresponded to a very noisy image with 
strong aliasing artifact. In contrast, under the rotating scheme of 
visiting 32 positions, the max g-factor (the fourth row) of 4-ele- 
ment RRFCA decreased to 1.6, which was comparable to the sta- 
tionary coil array with twice as many elements (max g= 1.7) in 
the second row. The reconstructed image was better with lower 
RMSE and artifact power [45] (RMSE=0.023, artifact 
power = 0.0082) compared to the 4-element stationary array 
(RMSE = 0.039, artifact power = 0.0231). However, fewer sensitiv- 
ity profiles were available when visiting fewer angular positions, 
leading to higher g-factors and less capability in imaging accelera- 
tion. As shown in the third row, the maximum g-factor increased to 
2.5 by only visiting six positions, and aliasing artifacts started to 
emerge in the reconstructed image. It is well known that the SNR 
is proportional to the square root of the number of channels 
[4,46]. To provide a fair comparison, an 8-channel RRFCA was sim- 
ulated by combining two sets of sensitivity profiles of the 4-chan- 
nel RRFCA with a 45° rotation. The g-map, reconstructed images 
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Fig. 7. Intra-subject case: comparisons between experimentally measured and numerically estimated sensitivity maps at position 2 (p2 in (a)) and position 3 (p3 in (b)). First 
row, experimentally acquired magnitude maps at p2; second row, numerically estimated magnitude maps at p2; third row: experimentally acquired phase maps at p3; fourth 


row, numerically estimated phase maps at p3. 


and the SNR maps were shown in the fifth row for such a coil array. 
Compared to the second row, we can see that with the same num- 
ber of coils, not only was the relative-SNR of RRFCA (relative- 
SNR = 2.83) higher than that of the stationary array (relative- 
SNR = 2.35), but also the global SNR map was more uniform. 


4. Discussion 
4.1. The library data 


The library data was used to provide sensitivity maps at all 
angular degrees, and the registration algorithm was used to bridge 
the discrepancy between the library and actual sensitivity maps. 
Theoretically, the library data can be either experimentally mea- 
sured or numerically calculated. The experimentally acquired sen- 


sitivity maps may be more realistic, but obtaining 3-dimensional 
maps at multiple positions can be time-consuming. Acquisition 
speed may be further restricted due to the potential heating prob- 
lem for ultra high field MRI. 

Commercially available electromagnetic (EM) software can be 
employed to calculate sensitivity maps and avoid such problems. 
With different algorithms, EM field distributions and related coil 
sensitivity can be calculated by solving Maxwell’s equations. The 
Method of Moments (MOM) [13,18,47-50], used for calculating 
the EM field, is efficient for homogeneous loads, but it is not feasi- 
ble to calculate heterogeneous dielectric loads due to the complex- 
ity of calculating the Green function [51]. The Finite Element 
Method (FEM) [52], discretising heterogeneous subjects into tetra- 
hedral or hexahedral elements, is capable of representing complex 
heterogeneous subjects smoothly and providing a more accurate 
solution. However, the FEM requires very large computational 
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Fig. 8. Inter-subject case: comparisons between experimentally measured and numerically estimated sensitivity maps at position 2 (p2 in (a)) and position 3 (p3 in (b)). First 
row, experimentally acquired magnitude maps at p2; second row, numerically estimated magnitude maps at p2; third row: experimentally acquired phase maps at p3; fourth 


row, numerically estimated phase maps at p3. 


resources for discretising the subjects and would have taken sub- 
stantial amount of time to accurately evaluate B4. The Finite-Differ- 
ence Time-Domain (FDTD) method [53-59] has advantages over 
the FEM as it simplifies the discretisation into regular boxes, and 
the iterative solution saves on computational resources. Working 
in conjunction with a graphical processing unit (GPU), the calcula- 
tion time can be dramatically reduced [60]. 


4.2. Transmit By 


A signal intensity image is determined by both transmit profile 
By and receive sensitivity By as illustrated in Eq. (7). At low field, 
the |By| is uniform and thus can be excluded from the sensitivity 
encoding matrix in the SENSE reconstruction. However, to include 
the |B} | influence in the SENSE reconstruction at ultra high fields, 


the actual sensitivity map in the encoding matrix should be in 
the form of Eq. (8) [14,61]. 


= Mo|B,"| sin (Vt|B; |) (7) 
sensitivityacuar = SİN (V)T|BY |snim) By (8) 


where y is the gyromagnetic ratio and Mo is proportional to the pro- 
ton density distribution. t and V denotes the RF pulse duration and 
coil driving voltage. The asterisk denotes the complex conjugate 
operation. 

In this work, an 8-element coil array was manufactured and 
actively detuned to provide an unchanged Bý in the course of rota- 
tion. However, the |By| inevitably changed as the array rotated, 
which may have introduced a small bias into the reconstruction 
when this changing field was not considered in the encoding 
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Fig. 9. g-map, image reconstruction and relative-SNR map comparisons between 4- and 8-element stationary coil array (first and second rows), 4-element RRFCA visiting 6 
positions (third rows), 4- and 8-element RRFCA visiting 32 positions (fourth and fifth rows). Left column, g-map comparisons; middle column, image reconstruction 
comparisons; right column, relative-SNR map comparisons. All SNR calculations were relative to that of a single rotating coil without undersampling (SNR = 1). 


matrix. To further improve the sensitivity encoding ability of the 
RRFCA and faithfully reconstruct the image, it may be advanta- 
geous to adopt a composite sensitivity concept [14]. Our future 
work will involve the development of the sensitivity estimation 
algorithm for the RRFCA working in transceive mode. 


4.3. SNR calculation 


As an important metric to quantitatively evaluate image qual- 
ity, the SNR calculation for MRI has been extensively studied 
{4,62-65]. At ultra high field (>3T), the SNR can be expressed as 
[62,66]: 


SNR Vsignal x B$ fo Wsin (Vyt|By |) |B; |dv (9) 
Vnoise V/Psample 


where the integration is performed over the volume of interest 
(VOI) and W is a weighting factor related to tissue and sequence. 


The Psampie is the power dissipated in the sample, which is too com- 
plex to be efficiently calculated for rotating array coils at ultra high 
fields. Thus the relative SNR was reported. In this work, to simplify 
the comparison, the Psample for different arrays were considered the 
same [67], since transmission power remained approximately the 
same for different arrays. 

In parallel imaging, the reduced k-space data are recovered by 
employing the coil sensitivity profiles. However, the noise is inev- 
itably amplified in the reconstruction process, especially with a 
high reduction factor. The SNR calculation when employing paral- 
lel imaging algorithms, such as SENSE, can be calculated as below 
[7], provided that the channel number does not change: 


SNR» 1 
SNRan gVR 
where g and R denotes the g-map and reduction factor, respectively. 


The SNR is also proportional to the square root of the number of 
channels [4,46]. In order to provide a fair SNR comparison, we 


(10) 
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simulated the 8-channel RRFCA by combining data obtained with a 
45° separation. 


4.4. RRFCA structure and data sampling 


4.4.1. Coil geometry 

Using a large number of RF coils can increase the SNR and sig- 
nificantly accelerate the imaging process [68-71]. However, plac- 
ing a large number of coils in a constrained space will decrease 
the coil size, leading to a shallower B, penetration. Namely, smaller 
coils receive less signal from the centre compared to larger coils. In 
addition, a coil array with higher density can increase the difficulty 
of decoupling. Coil coupling undermines the parallel imaging per- 
formance and reduces the SNR. By employing the rotating scheme, 
the RRFCA not only provides a large number of sensitivity profiles 
without additional RF channels, but also allows bigger coils to be 
used with less decoupling complexity. 


4.4.2. Data sampling 

The RRFCA is designed to sample one k-space line at each step- 
ping position, and then move to the next angular position for the 
following sampling. In this preliminary in vivo study, the stepping 
angle was manually adjusted and a full k-space matrix was sam- 
pled at each angular position. The accelerated acquisition was 
numerically modelled by extracting and combining the corre- 
sponding k-space lines. 


4.4.3. Rotating speed and acoustic noise 

In this proof-of-concept work, the RRFCA was manually rotated 
to sample data at several angular positions in order to validate the 
proposed in vivo sensitivity estimation algorithm. In our previous 
works [10,11], the rotating coil was pneumatically driven to 
achieve up to 870rpm for human head imaging with negligible 
acoustic noise, compared with noise generated by the gradient sys- 
tem. For animal imaging, the rotating speed can easily exceed 
10,000 rpm. Our recent experiments on a 9.4T pre-clinical system 
have shown that the rotating coil can faithfully reconstruct the 
images at speeds of, or exceeding, 5500 rpm. 


5. Conclusion 


In this paper, an in vivo rotational sensitivity estimation algo- 
rithm for the RRFCA was proposed and verified using human head 
imaging. By taking advantage of registration techniques and a 
library data, the algorithm was able to estimate sensitivity maps 
with arbitrary angular displacement. The estimated sensitivity 
maps were then fed into the rotating-SENSE algorithm to improve 
the imaging acceleration ability of the 4-element RRFCA. The 4- 
channel RRFCA outperformed the 4-element stationary array and 
was comparable to 8-element stationary PACs in terms of g-map 
and reconstruction quality. The 8-element RRFCA has been shown 
to significantly improve image quality compared to a stationary 8- 
element coil array. In the future, a sensitivity estimation algorithm 
for a transceive RRFCA with “composite sensitivity” will be devel- 
oped. By taking By into account in sensitivity encoding, the imag- 
ing acceleration ability of the RRFCA can be further improved. 
Additionally, an automatic rotation control system for the RRFCA 
will be developed by means of, for example, a non-magnetic 
piezo-electric or ceramic motor. 
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